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We develop an analytical approach for the description of an atomic Bose gas at unitarity. By 
focusing in first instance on the evaluation of the single-particle density matrix, we derive several 
universal properties of the unitary Bose gas, such as the chemical potential, the contact, the speed 
of sound, the condensate density and the effective interatomic interaction. The theory is also 
generalized to describe Bose gases with a finite scattering length and then reduces to the Bogoliubov 
theory in the weak-coupling limit. 
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Introduction. — Although strongly interacting sys- 
tems are known to be notoriously difficult to describe 
from first principles, remarkable progress in our under- 
standing has been made in recent years. For example, 
the universal nature of fermionic many-body systems 
with resonant two-body interactions has been success- 
fully studied experimentally [1]. This has been achieved 
by utilizing the high degree of control available in ul- 
tracold atomic gas experiments, which allow the in- 
vestigation of many-body systems from the weakly to 
the strongly interacting case, by using a magnetic-field- 
tunable Feshbach resonance [1-3] . The remarkable prop- 
erty of such resonant systems, which have an infinite scat- 
tering length and are therefore said to be at unitarity, is 
that at zero temperature there is no other length scale 
than the average interatomic distance that is set by the 
particle density n. As a result all thermodynamic quanti- 
ties, when appropriately scaled, can then be expressed in 
terms of a set of universal numbers. One of the most cru- 
cial quantities of the Fermi gas at unitarity is the chem- 
ical potential 

M=(l + /3)ei., (1) 

which is given by an universal constant times the Fermi 
energy ep = fi^k\/2ra, where kp — (67r^n/2s -I- 1)^/'^ is 
the Fermi momentum and s = 1/2 due to the hyperfine 
degrees of freedom. The universal constant /3 can be 
interpreted as describing the deviation from the ideal gas 
result due to interactions and was found to be /3 ~ —0.63 
experimentally as well as theoretically [4, 5]. 

Not only fermionic many-body systems are studied 
with ultracold atomic gases, but also its bosonic coun- 
terpart can be experimentally realized. In recent years, 
there has therefore been an increasing interest in the ex- 
perimental study of the strongly interacting Bose gas at 
low temperatures. It is expected on dimensional grounds 
that the Bose gas at unitarity, if stable, has similar uni- 
versal properties as that of the unitary Fermi gas. For 
instance Eq. (1) is expected to hold also but with s — Q 
and a different value of /3 due to the different statistics 
of the atoms. In contrast to fermionic cold atomic gases, 
the experimental study of bosons with increasing inter- 




Figure 1. The universal one-particle density matrix n(x)/n 
as a function fcF|x|, where the thin line is the contribution 
of the condensate and its phase fluctuations as mentioned in 
the text. From the difference of the graphs the contribution 
coming from all other fluctuations can be inferred. The con- 
densate density Uc/n is indicated on the right. 



action strength is complicated by the loss of atoms, as a 
consequence of a strong increase in the rate of inelastic 
three-body recombination processes caused by the ab- 
sence of the Pauli principle and the existence of Efimov 
trimers. As a result, up to now only a lower bound of 
/3 > —0.56 was determined experimentally [6]. These 
inelastic three-body processes result in the formation of 
molecules, which shows that the actual ground state of 
these gases is a Bose-Einstein condensate of molecules. 
Nevertheless, it may still be experimentally possible to 
create the meta-stable state of a Bose-Einstein conden- 
sate of atoms at large scattering lengths for a sufficiently 
long time. In view of this possibility, there has been con- 
siderable theoretical interest in the unitary Bose gas, but 
recent theoretical results strongly vary from /? ~ —0.34 
using renormalization-group techniques [7], /3 ~ 1.93 and 
P ~ —0.2 from variational (Jastrow) analyses [8, 9], to 
the complete instability of the Bose gas at unitarity using 
a Nozieres-Schmitt-Rink-like approach [10]. This clearly 
indicates that a systematic approach for the meta-stable 
atomic Bose gas at unitarity is still lacking. 
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In this Letter we present an analytical approach to the 
Bose gas at unitarity that can be improved systemati- 
cally. It is based on the realization that fluctuations in 
the phase of the Bose-Einstein condensate dominate the 
long-wavelength behavior of the system [11-13]. Going 
beyond the Bogoliubov theory, which is required at uni- 
tarity, these gapless phase fluctuations lead to infrared 
divergencies, which make it increasingly difficult to apply 
resummation procedures to find, for example, the effec- 
tive interaction at unitarity. Our approach circumvents 
these troublesome divergencies by exactly incorporating 
the phase fiuctuations, which will be confirmed by repro- 
ducing the exact form of the single-particle propagator in 
the long-wavelength limit as derived by Nepomnyashchii 
and Nepomnyashchii [14, 15]. The approach thus con- 
sists of isolating the phase fluctuations of the condensate, 
which is reminiscent of bosonization for fermions. In ad- 
dition, the theory is first renormalized due to all other 
fluctuations using the renormalization group. One of the 
outcomes of our theory, the universal one-particle den- 
sity matrix, is shown in Fig. 1. The thin line indicates 
the contributions to the density matrix from the conden- 
sate and its phase fluctuations, from which it is seen that 
these fluctuations are vital to the correct description of 
the Bose gas at unitarity. 

Renormalized Bosonization. — The Euclidean ac- 
tion of a Bose gas with a point interaction is given by 
S [(f)* , (/)] = / drdxC, where the Lagrangian density is 



riQ exp [lb 



2m 



1. 



(2) 



(/)(x,t) is the bosonic fleld, and the strength of the in- 
teraction is characterized by T^^ — Airafi^ /m with a the 
s-wave scattering length. To describe the Bose-Einstein 
condensate at zero temperature, we expand the field as 



(x, r) = y/uQ exp [i9{x, r)] -|- 0'(x, r), 



(3) 



where no should be viewed as the quasicondensate den- 
sity [3] and is not the density of atoms in the condensate 
tIq , as illustrated by Fig. 1. Roughly speaking, the first 
term of the expansion describes the low-energy modes 
of the field, as shown in Fig. 2, and includes the phase 
fiuctuations. The field 0'(x,r) describes the high-energy 
modes and is defined such that it does not contain phase 
fiuctuations and is thus orthogonal to the first term in 
Eq. (3). By inserting the expansion into Eq. (2), the ac- 
tion S [no, 9, 0'*, 0'] is obtained. There are several prop- 
erties of this action worth mentioning. 

First of all, the dominant low-energy physics is known 
to be due to the phase fiuctuations. It can be shown that 
the exact form of the phase- fiuctuation action S[9] in the 
long-wavelength limit is 



'(k,a;„) 



mc^/n 



{TiWnY + 2mc^ek; 



^(k,w„), (4) 
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Figure 2. Schematic representation of the expansion of the 
field (f) in terms of the condensate and its phase fluctuations 
and the non-phase fluctuations (j) , c.f. Eq. (3). 



where ujn = 2mT/h/3 is the Matsubara frequency, (3^^ — 
ksT is the temperature, k is the wave vector, et = 
/i^k^/2m is the atomic dispersion, and c is the speed of 
sound. In other words, Eq. (4) can formally be obtained 
from the full action S[na,9, cf)'* , cf)'] by integrating out all 
non- phase fiuctuations 4>' and the fiuctuations in tiq. 

Second, the accuracy of the action S[nQ,9, (j)'* .(j)'] can 
be improved systematically by incorporating the (j)' fluc- 
tuations into a renormalization of the action. However, it 
turns out to be more convenient to carry out this renor- 
malization at the level of S [(j)* , (j)] , Eq. (2) , and then ap- 
ply the expansion of the fleld, as in Eq. (3). The ex- 
act renormalization-group flow equation for the action 
S\6\6] is 



G' 



1^^ 



Here S [(j)* , (j>; A] is the effective action obtained by in- 
tegrating out all non-phase fluctuations above the mo- 
mentum hA, G' is the matrix propagator of the non- 
phase fluctuations, S-mt is the non-gaussian part of the 
effective action, the trace is over space, imaginary time 
and Nambu space $(k,a;„) = [(/)'(k,a;„), ^'*(— k, — a;„)]"^, 
and = S{k — A). Although there are no small pa- 
rameters in the theory of unitary Bose gases, the renor- 
malization group can distinguish between the relevance 
of the various coupling constants based on their scal- 
ing dimension under renormalization. As the effective 
interaction evaluated at zero momentum and zero fre- 
quency is expected to be a crucial variable, which how- 
ever also induces a flow of the chemical potential, let us 
here restrict our attention to these parameters. The run- 
ning of the chemical potential and effective interaction 
are in general found to be given by Ad/i/dA = I3^{^,g) 
and Adf//dA = (3g{fi,g). By solving these equations the 
renormalized action S [4>* , (f)', A] is found. Then, after in- 
serting the expansion of the fleld, the renormalized ac- 
tion 5 [riQ, 0, 0'*, 0'; A] is obtained. This action deflnes 
the propagator of the non-phase fluctuations in terms of 
the effective interaction. 

Lastly, to actually perform the above-mentioned renor- 
malization, the propagator of the non-phase fluctu- 
ations (j)' is needed. It is obtained by realizing 
that these fluctuations are determined from the action 
S[no,9,(j)'* ,(j)';A] with a non-fluctuating phase. The 
quadratic part of the action is then written as S [(j)'* , 0'] = 
5Ek.„*^(k,^n) [-^G'''(k,w„)] *(k,w„), where the 
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components of the symmetric 2x2 Green's function of 
the non-phase fluctuations are given by 



^,1 s fc if^n + eic + nog 



(k,a;„), 



-nog 



+ no (0r)(k,c^„), 



with the BogoHubov dispersion (/iWk)^ = ek(ek + 2no(7). 

The first term in the propagator is the standard BogoH- 
ubov propagator, which also includes the effect of phase 
fluctuations. The second term precisely removes these 
effects, which can for instance be identified by their pro- 
portionality to uq in the numerator. For consistency, 
we define mc^ = nog in accordance with the Bogoliubov 
dispersion. Notice, however, that the residues of the ex- 
act propagator of the phase fluctuations incorporate the 
renormalization of the superfluid density from no to n, 
which is not present in the Bogoliubov propagator. 

After the subtraction of the phase fluctuations, the 
propagator of the non-phase fluctuations is given by 



ihUn + Ck 



(5) 



and the anomalous averages vanish, i.e., {(t>'<j)') = 
= 0. The vanishing of the anomalous averages 
means that the Green's function G' is diagonal and this 
greatly simplifies the renormalization procedure of the 
interaction, which justifies our preference for solving the 
renormalization-group equation for S [0* , 0] . 

Before we turn to the determination of the effective in- 
teraction, we first show that our approach reproduces the 
exact propagator in the static long-wavelength limit de- 
rived by Nepomnyashchii and Nepomnyashchii, as men- 
tioned in the introduction. This is achieved using the 
one-particle correlation function ((/)(x, t)(/)*(0, 0)) , which 
equals no (exp [i (6i(x, r) - 61(0, 0))]) -I- (0'(x, t)0'*(O, 0)), 
and taking the Fourier transform. By expanding the ex- 
ponential, we find that the dominant long- wave length be- 
havior is due only to the first two terms in the expansion, 
where the first term is simply the exact phase-fluctuation 
propagator in Eq. (4) while the second term comes from 
a convolution of two propagators, namely 



h ^ (0(k,i:j„)(/i*(k,a;„)) 



-mc 



32^/2: 



3/2 



■log 



"-c/5V(5k,o'5„,o, 



which is the form of the exact propagator [14, 15] with Uc 
the condensate density defined later in Eq. (6) and V is 
the volume. The exact anomalous propagator is also re- 
produced, and is given by no (exp [i {9{x, r) -|- 9{0, 0))]). 
In particular, this leads to the counter-intuitive conclu- 
sion that the anomalous self-energy vanishes at zero mo- 
mentum and frequency [15]. 



Expand (l> Non-9uc. 9 
S [</>', <l>] *- S[na.e. (!>",(!>'] > S [no, 0, (4'*, <)>'] 

I ^ 1 

Effective interaction g and cliemical potential ^ due to fluctuations 

Figure 3. Schematic representation of the renormalization 
procedure that shows the self-consistent nature of our theory, 
as described in the text. 



To summarize and as illustrated in Fig. 3, the ac- 
tion S [(j)* , (j)] of the Bose gas can be systematically 
renormalized by the non-phase fluctuations (j)' using the 
renormalization-group flow equation, giving in particular 
rise to an effective coupling g and a renormalized chem- 
ical potential /i. The propagators of the non-phase fluc- 
tuations are determined self-consistently after expansion 
of the field. The theory includes the exact propagator of 
phase fiuctuations and can reproduce the exact propaga- 
tor {4>(f>*) in the long-wavelength limit. 

Universal Results. — Now that our theoretical frame- 
work has been established it is possible to determine, for 
example, the condensate density, the effective interaction 
and the chemical potential at unitarity. The conden- 
sate density can be calculated from the off-diagonal long- 



range order, namely Uc = lim 



|x|- 



(0(x,o)r(o,o)), 



no exp 



3/2 



2V2-7r 



nog 



||(67r2no)V3 



, (6) 



where it must be noted that an ultra-violet subtraction 
is needed in order to calculate the expectation value of 
the phase fluctuations using Eq. (4) with mc^ = not; 
[3] . This ultra-violet subtraction removes the ultra-violet 
divergences associated with a point interaction, and is a 
result of the renormalization of the bare coupling to T^^ . 
In order to determine the condensate density at unitarity, 
the quasicondensate density no needs to be eliminated in 
favor of the total density n = (0(x, t)0*(x, r)) using 



no 



(8^/2 -37r) f2m 



247r2 



[nog] 



3/2 



(7) 



where the second term is the contribution from the high- 
energy fluctuations n' = ((/)'(x, r)(/)'*(x, r)), see Eq. (5). 
As required, exactly the same ultra-violet subtraction 
was used for the high-energy fluctuations as in Eq. (6). 

To proceed, we must now determine the effective inter- 
action g in some approximation. Taking only the renor- 
malization of the coupling constant and the chemical po- 
tential into account, which turns out to be very accurate 
for the unitary Fermi gas [16], the beta functions are 
given by 



4 



47rA3 



2 471- A3 
(27r)3 



I"aI 



I^A 



2eA 

4 



(8) 



va\ 



1 

2el 



where the Bogohubov dispersion huj^ and the coherence 
factors |ukP = |vkP + 1 = (^(^k + Ck) /2/iWk are evalu- 
ated at A. The effective interaction is obtained by inte- 
grating its differential equation using the boundary con- 
dition g{A — oo) = T^^, where it must be noted that 
the effective interaction inside the Bogoliubov dispersion 
is the fully renormalized value ^(A — 0) which, as pre- 
viously explained, is determined self-consistently. Ulti- 
mately, we obtain 



f 
9 



1 



1 



4\/27r2 



2m 



3/2 



(9) 



Note that this equation can also be obtained directly 
as the result of a resummation of an infinite number of 
the diagrams shown in Fig. 4. The equation also shows 
the highly non-perturbative nature of the renormaliza- 
tion group. 

In the unitarity limit, T^^ — > oo, the effective inter- 
action and the condensate density in terms of the total 
density are found by solving Eqs. (6-9) to be 



ng ^ 

ep 32/3 

n 



(1 + A) 



1/3 



f.09, 



n 
n 



A 



1 + A 



0.31, 



/2V2-7r 
1 + A°^H^7^ 



0.55, 



where A = n' /uq = (8^2 - 37r) /3V2 ~ 0.45. The de- 
pletion from the condensate is given by 1 — ric/n ~ 0.45, 
which clearly differs from the density of particles con- 
tributing to the non-phase fluctuating modes n' by phase- 
fluctuation contributions. 

The universal one-particle density matrix n(x) = 
((^(x, 0)(/)*(0, 0)) is shown in Fig. 1 and can be seen to 
reduce in the large-distance limit to the condensate den- 
sity. Also the contribution due to the phase fluctuations 
is indicated, which terminates at ng/n at equal position. 
Another interesting property at unitarity is called the 
contact C and is related to the short-wavelength behavior 
of the one-particle density matrix, namely n(k) ~ C/k* 





Figure 4. The Feynman diagrams for the fluctuations 4>' that 
contribute to the beta function Pg and can be viewed as the 
building blocks for the ladder and bubble sums included in 
the Bethe-Salpeter equation for the effective interaction g. 




1/kpa 

Figure 5. The chemical potential as a function of scatter- 
ing length. The gray dashed line is the Bogoliubov chemical 
potential [3]. 



[17, 18]. The contact is 



C ^ fmgy ^ 1 1 

[2eFj 34/3(i_^;^)4/3 



0.14. 



The change in the chemical potential follows from inte- 
grating Eq. (8) and is given by A/i = 2n'g. According 
to the exact Hugenholtz-Pines theorem the chemical po- 
tential in our theory is then given by ^ = nog + Afi, such 
that the universal chemical potential is 



npg + 2n'g 



2 1 

327^ (iT 



I-2A 



1.42, 



which results in the universal constant /3 ~ 0.42. Fur- 
thermore, the speed of sound at unitarity is given by 



mc 



nog 



1 M 
1 + 2Aei^ 



0.53 



e-F 



0.75. 



The expected value for the speed of sound at unitarity 
in terms of the chemical potential is mc^ = n(d/i/dn) = 
2/i/3 ~ 0.66/x, which is remarkably close to our result and 
shows the accuracy of the simplest first approximation 
that we have presented here. 

Conclusions. — We have constructed an approach to 
describe Bose gases at unitarity which can be improved 
systematically by renormalization-group methods. As an 
outlook, all quantities can be found not only at unitarity, 
but it is possible using Eq. (9) to find all these quanti- 
ties from weak to infinitely strong coupling, as is shown 
for the chemical potential in Fig. 5. The generalization 
of the theory to frequency-dependent interactions and 
also non-zero temperature is straightforward and will be 
elaborated on in a future publication. Furthermore, we 
expect that the approach can be applied to other sys- 
tems with a broken continuous symmetry, where similar 
infrared divergencies occur as a consequence of the Gold- 
stone modes. We hope that our results stimulate fur- 
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thcr experimental developments toward unitarity-limited 
Bose gases in the near future. 
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